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Abstract 

The determination of \V u f,\ from inclusive semileptonic B decays is limited by uncertainties in mod- 
eling the decay distributions in b — > ulv transitions. The largest uncertainties arise from the limited 
knowledge of the appropriate b quark mass and Fermi momentum to use in the parameterization of 
the shape function. This paper presents a new method in which these shape function parameters 
are constrained by the same data used to measure | V u b \ . The method requires measurements of the 
momenta of both the charged lepton and the neutrino in semileptonic B decays. From these quan- 
tities two complementary observables can be constructed, one for discriminating between b — ► u£v 
transitions and background and the other for constraining the shape function. Using this technique 
the uncertaintites in \V u b\ from the shape function may be signficantly reduced. 
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1 Introduction 



The precise determination of \V u f,\ is of comparable importance to the measurement of CP asym- 
metries in probing the CKM sector of the Standard Model. The most precise determination of 
\V u b\ at present comes from a measurement of the endpoint of the charged lepton energy spectrum 
in inclusive semileptonic B decaysQ, which has a total relative uncertainty of 15%. The major 
uncertainties in the endpoint measurement arise from the background from b — > civ decays and the 
theoretical prediction of the fraction of the b — > ulv spectrum in the endpoint region. These two 
sources of error have opposite sensitivities to the lepton energy cut; raising the cut decreases the 
former and increases the latter. The total error is minimized with a lepton energy cut of ~ 2.2 GeV. 
Within the context of the lepton endpoint method, further progress requires an improved under- 
standing of the shape function that describes the b quark mass and momentum distribution in the 
B meson. 

A number of proposals have been made to reduce the theoretical uncertainties in the deter- 
mination of \Vyh\. These include measuring the recoil mass of the hadronic system in a b — ► ulv 
transition IH, measuring the invariant mass squared of the lepton pairQ or using a combination of 
variables^. In addition, the use of the photon energy spectrum in b — > 57 decays to reduce shape 
function-related uncertainties in the extraction of \V u b\ from semileptonic decays was proposed in 
Ref. m and applied in Ref. Q; more recent work|| suggests that the non- universality of the shape 
function may significantly limit the precision achievable with this approach. 

In this paper we discuss a technique by which the theoretical uncertainty in extracting \V u b\ 
from inclusive semileptonic decays at the T(4S) may be significantly reduced. The method requires 
the reconstruction of the neutrino momentum vector in addition to that of the charged lepton in 
semileptonic B decays. Q Using this information two quantities are measured, one directly sensitive 
to the shape function parameters and another offering good discrimination between b — * ulv and 
b — > civ decays. This approach using two complementary observables could be used in measure- 
ments of \V u b\ at present e + e~ B factories. 

2 The b — ► u£v Generator 

A Monte Carlo generator has been implemented to study b — > ulv decays. The generator is based 
on the triple differential decay width of Ref. ||, which is calculated including 0(a s ) corrections. 
The parton-level distributions are convolved with a shape function to obtain distributions in the 
experimental observables. The shape function can be written in terms of a single variable k + . The b 
quark mass appearing in the parton-level distributions is replaced by + k + . For the distribution 
of k + we take [|j 

F(k+) = N(l - K) a e {l+a)K ; k = ^± <1. (1) 

mB - m b 

The shape function can be parameterized using and a or, equivalently, using A = tub — rrib and 
Ai = -3A 2 /(1 +0). 



1 Technical issues involving neutrino reconstruction will not be discussed here; such reconstruction has already 
been used in analyses of exclusive semileptonic decays M. 
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3 Observables 



Detection of a charged lepton at an momentum high enough to allow identification is the starting 
point for any analysis of inclusive semileptonic decays. The T(4S I ) decays into BB, giving rise to 

B mesons that are nearly at rest, with velocity j3 = ^Jl — (2mB/?7ir(4S)) 2 — 0.06. The lepton 
momenta in the B meson rest frame are smeared by the B motion in the T(4S) rest frame. The 
invariant masses of the lepton pair and of the recoiling hadronic system are unaffected by the boost. 

Two types of observables are constructed. The first type allows discrimination between b — * u£v 
decays and background. For this purpose, the invariant mass of the hadronic recoil, m^, can be 
used, provided the other B meson in the event is fully reconstructed. Alternatively, a measurement 
of the missing momentum, which can be obtained without fully reconstructing the other B meson, 
can be combined with the charged lepton momentum to determine q 2 , the invariant mass squared 
of the lepton pair. The measurement of q 2 , in contrast to the neutrino energy itself, does not suffer 
from the unknown direction of the parent B meson in the T(4S) frame. Using the charged lepton 
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Figure 1: q 2 vs. Eg in the T(4S) rest frame. The diagonal line and the curve are contours for 
s max _ q an( ^ m 2^ respectively. The left plot shows b — ► civ transitions and the right plot shows 
b — > xx.lv transitions. 



energy and q 2 one can define a variable to discriminate between b — > xxlv and b — ► civ decays. In 
the B meson rest frame, ignoring lepton masses, the neutrino energy satsifies E u > q 2 /4Ei. Setting 
the neutrino energy equal to its minimum value determines the maximum kinematically allowed 
invariant mass squared of the hadronic recoil system: 

s ™ ax = ml + q 2 - 2m B E e - 2m B (B rest frame) (2) 
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Accounting for the boost into the T(4S I ) frame and respecting the condition m/j + < m-Q gives 



2 i 2 

m B +q 



2m.BE i 



l±(3 




for ± Ei > ±i (m B - m re f) W , 



(3) 



ra B - m re f 



otherwise. 



In the above Eg is the charged lepton energy in the T(4S) rest frame. An appropriate choice of 
reference mass for the definition of s™ ax is Md. Events containing b — > u£u decays are selected by 
requiring s™ ax < m 2 ut , with m cut — m B ; no event with a true hadronic recoil mass above m cu t can 
enter the sample unless q 2 (or, less likely, Ei) is misreconstructed. Constant values of s™ ax map 
out contours in the q 2 -Ei plane as shown in Fig. [I]. Decays from the dominant b — > c£u process 
cannot contribute in the region s™ ax < m 2 ^ unless they are mismeasured. 

As seen in Fig. pL the shape of the b — > u£u spectrum in m\ is sensitive to both the assumed 
b quark mass and to the Fermi momentum. The spectrum in s™ ax depends on but is largely 
insensitive to the Fermi momentum; given that it is based largely on q 2 , this is as expected. In 
these and subsequent plots the charged lepton energy is required to exceed 1 GeV in the T(4S) 
frame to ensure that it can be cleanly identified.^ 
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Figure 2: The effect of varying the shape function parameters on the distributions of and s™ ax 
for b — > uiiy decays. In each case the charged lepton exceeds 1 GeV in the T(4S) frame. The 
b — > u£v signal region is to the left of the dashed vertical line. 



2 The efficiency of this requirement is 87% and is not very sensitive to the shape function. 
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The second type of observable is sensitive to the effective b quark mass. In the rest frame of 
the b quark, is related to the W properties: 



w ~~ 2mb ~ 

where m u is the effective mass of the final state u quark and the equality in the last relation holds in 
the limit m u <C my,. The quantity E^ + \pff\ is an estimator for the effective b quark mass. The 
boost from the b quark rest frame into the T(4S) rest frame smears the estimate but introduces 
little bias, since the boost direction is uncorrelated with the W direction and 7 ~ 1. The average 
value {Ew + |pw|), measured in the T(4S) rest frame from the lepton momenta, is an observable 
with substantial sensitivity to the effective b quark mass. 

The strategy for determining \V u b\ uses either or to separate b — > w.lv decays from 
background and (£\y + |pW|) to limit the variation in shape function parameters that must be 
considered. 

How well can (£\y + \pw\) be measured? The r.m.s. of the generated £\v + |pw| distribution 
is ~ 0.25 GeV after cuts are applied on s™ ax or to select b — > ulv events. The experimental 
resolution on missing energy will further broaden this; the r.m.s. missing energy resolution in 
the e + e~ B factories is 0(0.2 GeV). With B-factory data samples the statistical uncertainty in 
(£\y + \pw\) can be reduced below ~ 0.02 GeV. However, it will be an experimental challenge to 
understand systematic biases^] at this level. 

4 Extracting \V u b\ 

The two parameters and a of the shape function are varied to estimate the change in the fraction 
of b — ► xtlv decays in the signal region. Fig. || shows the change in the fraction f u of b — » u£u 
transitions in the signal region versus (£\y + |pw|) for different choices of shape function parameters. 
The variations m h = (4.8 ± 0.15) GeV and a = 1.29+J;|J correspond to A = (0.48 ± 0.15) GeV 
and Ai = (— 0.3^ ' 20 ) GeV 2 , respectively, and are a reasonable a priori choice at present. The 
observable (E\y + \pw\) has sensitivity to both and a; the sensitivity to m-b is reduced (i.e. 
A(£\y + |pw|)/Amb < 1) due to the cuts used to reject b — > civ decays. 

With no measurement of {Ey? + |pw I } the uncertainty on | V u \, \ due to the shape function (about 
~ 9% for the requirement s™ ax < 3.2 GeV 2 ) is dominated by the uncertainty in m^. A simultaneous 
measurement of {E\y + |pw|) could significantly reduce this uncertainty. For example, achieving 
a precision of 0.04 GeV on (£\y + \pw\) reduces the uncertainty on \V u b\ to ~ 5% for a cut on 
s max ^ 3.2 GeV 2 (see the shaded band in Fig. ||). Similar results are obtained if < 1.7 GeV is 
used to select b — > uiu decays. If a stiff cut on the lepton energy is required in addition to the 
cut on s h Iiax , or if the cut on is reduced to 1.5 GeV, the theoretical error is increased and the 
ability to reduce it with a given precision on (£\y + |pw|) is decreased. However, a measurement 
of (£\y + \pw\) would provide an important cross-check on the externally chosen parameter values 
in all cases. 



3 Biases not properly accounted for in Monte Carlo simulations. Note, however, that uncertainties on the missing 
energy resolution do not affect (Ew + \pw\) directly. 
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Figure 3: The ratio / u (m.b, a)//u(4.8, 1.29) is shown versus (I?w + \pw\) for different choices of the 
parameters nib and a, with the m\, values indicated by the numbers on the plot. The leftmost, center 
and rightmost curves correspond to choices a = 0.38, 1.29 and 3.60, respectively. The requirement 
s max 3.2 GeV 2 is used in (a); in (b) the additional requirement Ei > 2.1 GeV is made. In (c) and 
(d) the requirements < 1.7 GeV and < 1.5 GeV are made, respectively. In each case the 
charged lepton energy is required to exceed 1 GeV in the T(4S) frame. The shaded bands show an 
example of the impact a measurement of (£\y + |pw|) with 0.04 GeV precision would have. The 
values of / u (4.8, 1.29) corresponding to requirements (a)-(d) are, respectively, 0.22, 0.15, 0.67 and 
0.57. 
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5 Discussion 

The use of (E\y + Ww I ) to constrain the shape function is analogous to the use of the photon energy 
spectrum from b — > 57 decays and amounts to using the virtual W to probe the b quark decay. 

In order to use the b — > S7 spectrum to constrain the uncertainty on \ V u b\, the 7 spectrum must 
be measured to as low an energy as possible and deconvoluted to obtain an estimate of A (and, 
perhaps eventually, Ai). Care must be taken to ensure that the theoretical expressions used to 
extract A from b — > S7 and to determine the dependence of \V u b\ on A are compatible. Corrections 
to the universality of the shape function must also be considered when applying this information 
to semileptonic b — > u£u decays. 

There are several advantages to the approach described in this paper. Since the same event 
generator is used for extracting jV^I and determining the relationship between + |pW|) and 
the shape function, no deconvolution is necessary. The need for a common theoretical framework 
for extracting and using the shape function parameters is trivially met; the parameters of the shape 
function need not even be extracted. 

Quantifying the shape function parameters is of interest in its own right. In the absence of 
experimental cuts, the perturbative relationship between A and (£?w + IjPwI) 111 the B meson rest 
frame can be determined (to 0(a s )) from formulae in Ref. ||: 

A = m B - (E w + |pw|) - — - _ . (4) 

The relationship between A and the measured {Ey? + |pw|) needs to be calculated in the presence 
of the experimental cuts required to reject b — » civ events. Extracting Ai is not feasible due to the 
large intrinsic width of the distribution and direct sensitivity to the missing energy resolution. 

The method outlined in this paper holds promise for reducing the theoretical uncertainty in 
determining the fraction of b — > ulv decays in the experimentally accessible region. The s™ ax 
variable offers a way of extracting a clean signal for b —* u£v decays with modest theoretical 
uncertainties. While it has a lower acceptance for b — > uiv decays than does the recoil mass m^, it 
may be advantageous experimentally because it does not require full reconstruction of the second 
B meson in the event. The use of (£w + |Pw|) ma Y reduce the range of parameter space that 
need be considered in evaluating the uncertainty on \V u b\, and provides an important cross-check 
on externally motivated parameter choices. The ultimate experimental precision on \ V u b\ is hard to 
predict at present, but a precision approaching 5% can be envisaged. There is strong motivation 
to pursue such a measurement. 
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